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Abstract 

We compare the coior-evaporation model (CEM) and nonrelativistic QCD (NRQCD) factoriza- 
tion predictions for inclusive quarkonium production. Using the NRQCD factorization formulas 
for quarkonium production and for perturbative QQ production, we deduce relationships that are 
implied by the CEM between the nonperturbative NRQCD matrix elements that appear in the 
factorization formula for quarkonium production. These relationships are at odds with the phe- 
nomenological values of the matrix elements that have been extracted from the Tevatron data for 
charmonium production at large transverse momentum. A direct comparison of the CEM and 
NRQCD factorization predictions with the CDF charmonium production data reveals that the 
CEM fits to the data are generally unsatisfactory while the NRQCD factorization fits are gener- 
ally compatible with the data. The inclusion of Jzt smearing improves the CEM fits substantially, 
but significant incompatibilities remain. The NRQCD factorization fits to the Xc data indicate 
that multiple gluon radiation is an essential ingredient in obtaining the correct shape of the cross 
section as a function of pr- 
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I. INTRODUCTION 



In general, the production rates for specific hadronic states in high-energy processes are 
difficult to understand from first principles because they involve nonperturbative aspects of 
QCD in an essential way. The inclusive production rates for heavy-quarkonium states have 
two aspects that facilitate their understanding from first principles. First, the mass itlq of 
the heavy quark and the antiquark that are the constituents of the quarkonium are large 
compared to the scale Aqcd that is associated with the nonperturbative aspects of QCD. 
Second, the inclusive nature of the quarkonium production process may make it less sensitive, 
at transverse momenta pr that are much greater than Aqcd, to nonperturbative effects that 
are associated with the formation of color-singlet hadrons from the colored partons. 

Two models for the inclusive production of heavy quarkonium were introduced in the 
1970's: the color-singlet model (CSM) and the color-evaporation model (CEM). These mod- 
els make very different assumptions about the roles played by the colors and spins of the 
heavy quark and antiquark in the production process. In the CSM, it is assumed that a 
specific quarkonium state can be formed only if the QQ pair is created in a color-singlet 
state with the same angular-momentum quantum numbers as the quarkonium. In the CEM, 
the probability of forming a specific quarkonium state is assumed to be independent of the 
color of the QQ pair. In some versions of the CEM, the probability of forming a specific 
quarkonium state is also assumed to be independent of the spin of the QQ pair. In spite of 
the very different assumptions that are made in the CSM and CEM, both models enjoyed 
considerable phenomenological success through the 1980's and into the 1990's. 

In 1995, a new approach for describing inclusive heavy-quarkonium production based on 
first principles was developed: the nonrelativistic QCD (NRQCD) factorization approach 
[1]. It makes use of an effective field theory called NRQCD [2, 3] to exploit the large 
mass of the heavy quark, and it further assumes that, owing to the inclusive nature of 
the production process, traditional factorization methods can be exploited to establish a 
factorization formula. The NRQCD factorization approach incorporates aspects of both the 
CSM and CEM and can be regarded as a unification of these two models within a consistent 
theoretical framework. It can be summarized by the NRQCD factorization formula, which 
separates short-distance, perturbative effects involving momenta of order ijiq from long- 
distance, nonperturbative effects. 

The nonperturbative factors in the NRQCD factorization formula are NRQCD matrix 
elements. NRQCD can be used to give rough predictions of the relative sizes of the matrix 
elements. These predictions are based on the leading power behavior of the matrix elements 
as a function of the typical heavy-quark velocity v in the quarkonium rest frame. In order to 
obtain an estimate, one makes the assumption that, for each matrix element, the coefficient of 
the leading power of v is of order unity. This last assumption may not be reliable. However, 
one can ignore this f-scaling information and treat the NRQCD factorization formula as a 
general phenomenological framework for analyzing inclusive heavy-quarkonium production. 
Any model that can be described in terms of QCD processes at short distances, including the 
CSM and the CEM, can be formulated in terms of assumptions about the matrix elements 
in the NRQCD factorization formula. 

In this paper, we derive relationships between the NRQCD nonperturbative factors that 
follow from the model assumptions of the CEM. We find that these relationships are often 
poorly satisfied by phenomenological values of the NRQCD matrix elements. Furthermore, 
the relationships sometimes violate the v -scaling rules of NRQCD. We conclude that the 
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CEM and NRQCD provide very different pictures of the evolution of a heavy quark- ant iquark 
pair into a quarkonium. 

Using existing data, one can exclude the CSM as a quantitative model of heavy- 
quarkonium production. In 1995, the CDF collaboration measured the cross sections for 
the prompt production of J/ip and ip{2S) in pp collisions at a center-of-mass energy of 
1.8 TeV. It discovered that the cross sections are more than an order-of-magnitude larger 
than those predicted by the CSM. This dramatic discrepancy eliminated the CSM as a 
viable model for inclusive heavy quarkonium production. (A detailed discussion of these 
comparisons between the CSM and the experimental data and references to the published 
experimental results and theoretical results can be found in Ref. [4].) 

The CEM can also be ruled out on the basis of the simple qualitative prediction that the 
ratio of the inclusive production rates for any two quarkonium states should be independent 
of the process. The most dramatic violation of this prediction that has been observed is in the 
fraction of J/^'s that come from decays of the P-wave charmonium states Xci an d Xc2- This 
fraction is measured to be 0.11 ±0.02 in B decays and 0.297±0.017 stat ±0.057 syst for prompt 
production at the Fermilab Tevatron. The version of the CEM in which the probability for 
the formation of a quarkonium is assumed to be independent of the spin state of the QQ 
pair can be ruled out on the basis of several other simple qualitative predictions. One such 
prediction is that the inclusive production rate of a quarkonium state should be independent 
of its spin state, so that it should always be produced unpolarized. This prediction is 
contradicted by the observation of nonzero polarization of J/^'s in e + e~ annihilation at the 
B factories and by the observation of nonzero polarization of the bottomonium states T(2S) 
and T(3S) in a fixed-target experiment. A further prediction of this version of the CEM is 
that the production rates for the P-wave charmonium state Xcj should be proportional to 
2j + 1, and, hence, that the ratio of the inclusive cross sections for Xci and x& production 
should be 3/5. For prompt production at the Tevatron, this ratio has been measured to be 
1.04 ± 0.29 sta t ± 0.12 syst . (A detailed discussion of these comparisons between the CEM and 
the experimental data and references to the published experimental results and theoretical 
results can be found in Ref. [4].) 

Since the CEM is only a model, it can be salvaged simply by declaring it to have a limited 
domain of applicability. The failure of the predictions for polarization can be avoided by 
declaring the model to apply only to cross sections that are summed over the spin states 
of the quarkonium. In the case of the predictions that ratios of quarkonium cross sections 
should be the same for all processes, the most dramatic failures can be avoided by declaring 
the model to apply only when the total hadronic energy is sufficiently large. This condition 
can be used to exclude applications to B decays and to e + e~ annihilation at \fs =10.6 GeV. 
In this case, the CEM reduces essentially to a model for inclusive production of quarkonium, 
without regard to spin, in high-energy fixed-target experiments, pp collisions at the Tevatron, 
and pp collisions at the LHC. 

It may be possible to exclude even this limited version of the CEM, given sufficiently ac- 
curate experimental information. In this paper, we present a quantitative comparison of the 
predictions of the CEM and the NRQCD factorization approach for transverse-momentum 
distributions of heavy quarkonia at the Tevatron. We restrict our attention to the region of 
transverse momentum comparable to or larger than the quarkonium mass, where the effects 
of multiple soft-gluon emission are not so important. We find that the NRQCD factorization 
approach gives a significantly better fit to the available data. 

The remainder of this paper is organized as follows. In Sec. II, we describe the NRQCD 
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factorization formula for quarkonium production and the f-scaling rules for the NRQCD 
matrix elements. In Sec. Ill, we describe the CEM, use the NRQCD factorization formula 
to derive expressions for the ratios of NRQCD matrix elements that are implied by the CEM 
assumptions, and deduce the f-scaling rules that follow from these ratios. Sec. IV contains 
a comparison of the ratios implied by the CEM with the ratios of phenomenological matrix 
elements that have been extracted from the Tevatron data for J /if) and i[)(2S) production. 
Sec. IV also contains a direct comparison of fits of the CEM predictions and the NRQCD 
factorization predictions to the Tevatron data. Sec. V contains similar comparisons for Xc 
production. Finally, we present our conclusions in Sec. VI. 

II. NRQCD FACTORIZATION FORMULA 

The NRQCD factorization formula for the inclusive cross section for production of a 
specific heavy-quarkonium state H is 

a[AB ^H + X] = J2 c* B (A)(0»(A)). (1) 

n 

Here, A and B are light hadrons, photons, or leptons, and A is the ultraviolet cutoff of the ef- 
fective theory. The c^ B are short-distance coefficients that can be calculated in perturbation 
theory by matching amplitudes in NRQCD with those in full QCD. The matrix elements 
(O^) are vacuum-expectation values of four-fermion operators in NRQCD, evaluated in the 
rest frame of the quarkonium. These matrix elements contain all of the nonperturbative 
physics of the evolution of a QQ pair into a quarkonium state. The operators have the form 

O" = X ^ n ^V H (A)^K> nX , (2) 

where if) is the two-component (Pauli) spinor that annihilates a heavy quark, \ is the two- 
component spinor that creates a heavy antiquark, and V H is a projector onto states that in 
the asymptotic future contain the quarkonium H plus light partons X whose energies and 
momenta lie below the cutoff A of the effective field theory: 

V H (A) = ^2\H + X,t^oo)(H + X,t^oo\. (3) 
x 

The factors K n and n' n in Eq. (2) are direct products of a color matrix (either the unit matrix 
or the matrix T a with octet index a), a spin matrix (either the unit matrix or the matrix 
a 1 with triplet index i), and a polynomial in the QCD covariant derivative D = V + igA 
and the QCD field strengths. The NRQCD factorization formula in Eq. (1) was proposed 
in Ref. [1]. Some hard-scattering factorization formulas, such as those for deep-inelastic 
scattering, Drell-Yan lepton-pair production, and e + e~ annihilation into hadrons, have been 
proven to hold to all orders in the strong coupling a s . The derivation of the NRQCD 
factorization formula is not at this level of rigor. 1 However, in this regard it is no different 



A recent study of certain two-loop contributions to quarkonium production [5] has revealed that, if fac- 
torization is to hold, then the color-octet NRQCD production matrix elements must be modified from the 
form given in Eq. (2) by the inclusion of light-like cikonal lines that run from each of the QQ bilinears to 
the far future. It is not known if this modification preserves the factorized form in higher orders. 
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TABLE I: Velocity-suppression factors for NRQCD matrix elements in 5-wave and P-wave QQ 
channels in NRQCD and in the CEM. The 1 or 8 indicates the color channel and 2s+1 Lj indicates 
the angular-momentum channel. For NRQCD, the ^-suppression factors up to order v 4 are given 
for representative S'-wave and P-wave multiplets. For the CEM, the orders of the ^-suppression 
factors are independent of the quarkonium state H, as is described in Sec. III. 





1,^0 l, 3 5i S, 1 ^ 8, 3 S 1 


l/Pi 1, 3 P l, 3 Pi 1, 3 P 2 8/Pi 8, 3 P 8, 3 P! 8, 3 P 2 


NRQCD Factorization 


j/y 


1 V A V 3 

1 v 3 v A 


4 4 4 
V V V 


h c 

XcO 
Xcl 
Xc2 


v 2 

v 2 
v 2 
v 2 


v 2 

v 2 

v 2 

v 2 


Color-Evaporation Model 


H 


1111 


v 2 v 2 v 2 V 2 v 2 v 2 v 2 v 2 



from the factorization formula for semi-inclusive production of a hadron in hadron-hadron 
collisions. The existing all-orders proofs of factorization formulas require that the observed 
scattered particle be produced at a large transverse momentum compared with the QCD 
scale Aqcd and that the cross section be sufficiently inclusive. "Sufficiently inclusive" means 
that the variables in which the cross section is differential cannot assume values that restrict 
final-state parton momenta in the parton-level cross section to be within order Aqcd of soft 
or collinear singularities. 

The matrix elements in Eq. (1) fall into a hierarchy according to their scaling with the 
velocity v of the heavy quark (or antiquark) in the quarkonium rest frame, v 2 ~ 0.3 for 
charmonium, and v 2 ~ 0.1 for bottomonium. In practice, the summation over these matrix 
elements is usually truncated at a low order in v. The NRQCD factorization formalism 
has enjoyed a good deal of phenomenological success in describing inclusive quarkonium 
production at hadron, ep, and e + e~ colliders and in fixed-target experiments. 2 

A standard set of the NRQCD operators 0% that appear naturally in cross sections that 
are summed over the spin states of the quarkonium was introduced in Ref. [1]. They are 
denoted by ( 2s+1 Lj) and (!?| / ( 2s+1 P,), where the subscript indicates the color state of the 
QQ pair (1 for singlet and 8 for octet), and the argument indicates the angular-momentum 
state of the QQ pair (s is the total spin quantum number, L = S, P, . . . indicates the 
orbital-angular-momentum quantum number, and j is the total-angular-momentum quan- 
tum number). There are implied sums over the spin states of the quarkonium H. 

The velocity-scaling rules of NRQCD imply an intricate pattern of suppression factors 
for the NRQCD matrix elements (0%). The suppression factors, up to order t> 4 for S'-wave 
and P-wave QQ channels, are given in Table I for the representative S'-wave multiplet that 
consists of the charmonium states r] c and J / ip and for the representative P-wave charmonium 
multiplet that consists of the charmonium states h c , Xco, Xa, and Xc2- 



2 See Ref. [4] for a recent summary of the phenomenology of quarkonium production. 
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III. COLOR-EVAPORATION MODEL 



The color-evaporation model (CEM) was first proposed in 1977 [6-9]. In the CEM, the 
cross section for production of a quarkonium state H is some fraction F H of the cross section 
for producing QQ pairs with invariant mass below the MM threshold, where M is the lowest- 
mass meson containing the heavy quark Q. The fractions Ffj are assumed to be universal 
so that, once they are determined by data, they can be used to predict the cross sections in 
other processes and in other kinematic regions. The cross section for producing QQ pairs 
that is used in the CEM has an upper limit on the QQ pair mass, but no constraints on the 
color or spin of the final state. The QQ pair is assumed to neutralize its color by interaction 
with the collision-induced color field, that is, by "color evaporation." In some versions of the 
CEM [10], the color-neutralization process is also assumed to randomize the spins of the Q 
and the Q. The CEM parameter F H is the probability that a QQ pair with invariant mass 
less than 2mju, where tum is the mass of the meson M, will bind to form the quarkonium H. 
That probability is assumed to be if the QQ pair has invariant mass greater than 2m m- 

In the CEM, the production cross section for the quarkonium state H in the collisions of 
light hadrons, photons, or leptons A and B is 

r 4m li Ary 

& cem [AB ^H + X] = F H dm% Q 1 ~^[AB -> QQ + X], (4) 

J 4m 2 dm QQ 

where mqq is the invariant mass of the QQ pair, m is the heavy-quark mass, and da/dm^Q 
on the right side is the inclusive differential cross section for a QQ pair to be produced 
in a collision of A and B. There is an implied sum over the colors and spins of the final- 
state QQ pair. This is where the central model assumptions of color evaporation and spin 
randomization manifest themselves. 

If A and/or B are hadrons or photons, the cross section for AB — > QQ + X can be 
expressed as convolutions of parton distributions for A and/or B and a parton cross section. 
At leading order in a s , the parton process ij — > QQ creates a QQ pair with zero transverse 
momentum, and the differential cross section dajdp\ is proportional to 8{p\). At next-to- 
leading order in a s (NLO), there are parton processes ij — > QQ + k that create a QQ pair 
with nonzero pr- The complete NLO differential cross section is a distribution that includes 
singular terms proportional to 5{p"j) and but whose integral over p\ is well behaved. 

Some kind of smearing over p? is necessary to obtain a smooth p? distribution that can be 
compared with experiment. The physical origin of the smearing is multiple gluon emission 
from the initial- and final-state partons. A rigorous treatment of the effects of multiple gluon 
emission requires the resummation of logarithmic corrections to all orders in a s [11-13]. A 
simple phenomenological model for the effects of multiple gluon emission is kx smearing, 
in which the colliding partons are given Gaussian distributions in the intrinsic transverse 
momentum with a width (kj<) that is treated as a phenomenological parameter. 

Complete NLO calculations of quarkonium production in hadronic collisions using the 
CEM have been carried out in Refs. [14, 15], using the exclusive QQ production code of 
Ref. [16] to obtain the QQ pair distributions. There are also calculations in the CEM beyond 
LO that use only a subset of the NLO diagrams [10] and calculations that describe the soft 
color interaction within the framework of a Monte Carlo event generator [17]. Calculations 
beyond LO in the CEM have also been carried out for jp, 77, and neutrino-nucleon collisions 
and for Z decays [18-22]. 
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We now proceed to elucidate the relationship between the CEM and the NRQCD fac- 
torization formula. According to the NRQCD factorization formalism, the differential cross 
section for the process AB — > QQ + X is given by 

-^[AB^QQ + X] = J2c^J2J2 /^(e* fc VK«-4(fc 2 + m 2 )], (5) 

QQ n spins colors 

where the c^ B are the same short- distance coefficients that appear in Eq. (1). We have 
suppressed the dependence of the coefficients and the operators on the ultraviolet cutoff A 
of the effective field theory. The operator £)^( +fc ) < 2(~ fc ) j s analogous to the operator in Eq. (2), 
except that the quarkonium state H is replaced by a perturbative state that consists of a 
Q and a Q with momenta ±fe and with definite spin and color indices that have been 
suppressed. Those suppressed indices are summed over in Eq. (5). It is convenient to define 
a QQ operator that includes the sum over the colors and spins and an average over the 
directions of the momenta ±k of the Q and Q : 

opik) = (7 ^EE ^ (+fc)0( - fc) ) ^<x, (6) 

\ spins colors / 

where dflk is the element of angular integration of k and k = \k\. 

Comparing Eqs. (1) and (4) and making use of Eq. (5), we see that the CEM implies that 

__ 1 /*&max _ 

T, C n B (°n)=F H ^n B ^ jf ^dk(0^(k)), (7) 



where k max = \fm 2 M — m 2 . Eq. (7) is the central relation that connects the NRQCD factor- 
ization approach with the CEM. If NRQCD factorization can be established to all orders in 
perturbation theory, then Eq. (7), which is a statement about the perturbative production of 
QQ states, must hold rigorously. This is to be contrasted with the NRQCD factorization ap- 
proach itself, in which nonperturbative effects could conceivably spoil the factorization, even 
if the factorization formula can be established to all orders in perturbation theory. Note, 
however, that there is an implicit assumption in Eq. (7) that the sum over n converges. The 
sum need not converge rapidly in order to establish the relationship between the NRQCD 
factorization approach and the CEM, but it must converge. The expansion parameter in the 
sum is k^Jrr? ~ v 2 . Hence, if the sum is to converge, we must have k^Jm 2 < 1. In the 
charmonium system, M is the D meson, with mass 1.86 GeV. Taking m = 1.5 GeV, we find 
that ^ ax /m 2 w 0.54. In the bottomonium system, M is the B meson, with mass 5.28 GeV. 
Taking m = 4.7 GeV, we find that k^Jm 2 m 0.26. Thus the assumption k 2 m& Jm 2 < 1 is 
satisfied for the bottomonium system, but is only marginally satisfied for the charmonium 
system. There is also an implicit assumption that the energy of the produced quarkonium 
is not near the kinematic limit for the process. Near the kinematic limit, the v expansion of 
NRQCD breaks down, and, so, one needs to carry out a resummation of classes of NRQCD 
matrix elements in order to maintain the accuracy of the calculation [23] . Near the kinematic 
limit, the short- distance coefficients c£ B contain large logarithms that must be resummed 
as well [24]. 

Now we wish to establish that the equality in Eq. (7) must hold independently for each 
term in the sum over n. If we truncate the series at a finite number of terms, as is often done 
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in the phenomenology of quarkonium production, then the short-distance coefficients c„ 
generally vary independently of each other as the kinematic variables of the incoming and 
outgoing particles vary. However, in some processes, for example, in quarkonium production 
as a function of p T at the Tevatron, some of the short-distance coefficients show identical 
or nearly identical behavior as a function of the kinematic variable(s). Then one can estab- 
lish the equality only of linear combinations of matrix elements. However, one can make 
a stronger assumption that the CEM must hold, not just for a particular quarkonium pro- 
duction process, but for all possible quarkonium production processes. Then, since Eq. (7) 
must hold for an arbitrarily large number of processes that have different short-distance 
coefficients, the equality in Eq. (7) must hold independently for each term in the sum over 
n. Under this assumption, the CEM predicts that the NRQCD production matrix elements 
are given by 

{O " ) = 2^ Fh J k2dk (°n Q (k)). (8) 

The matrix elements on the right side of Eq. (8) can be computed in perturbation theory. 
Their dependences on k are governed by the powers of the covariant derivative D that 
appear in the factors K n and k' u of the operators. For the matrix element of leading order 
in v that corresponds to an operator of orbital-angular-momentum quantum number /, the 
matrix element is proportional to k 21 . Hence, the integral on the right side of Eq. (8) is 
proportional to k^£/(2l + 3). Since A; max scales as v, the matrix element is suppressed 
as v 21 compared to the matrix element of an S*-wave operator. Thus, the CEM implies 
a velocity-suppression pattern that is independent of the quarkonium state and depends 
only on the orbital-angular-momentum quantum number of the QQ pair. The suppression 
pattern for S'-wave and P-wave matrix elements is shown in Table I. It should be contrasted 
with the intricate suppression pattern implied by NRQCD. We note that the relation (8) 
partially satisfies the constraints that are imposed by the velocity-scaling rules of NRQCD, 
in that there is an additional power of v 2 for each unit of orbital angular momentum in the 
operator. However, powers of v that do not arise from the orbital angular momentum are 
not reproduced in the relation (8). This fact has been noted previously by Beneke [25]. 

Now let us state explicitly the relationships between the NRQCD matrix elements that 
are implied by the CEM relation (8). In the operator O n , let the subscript n represent the 
angular-momentum quantum numbers (s, /, and j) and the color state (singlet or octet). 
The standard forms of the operators are given in Ref. [1] for the first few S- and P-wave 
channels. The spin-triplet operators are normalized so that, if one makes the replacement 
<7j <E> (Tj — > Sij(l <8 1), they become (2j + 1) / (2Z + 1) times the corresponding spin-singlet 
operators. The color-octet operators are normalized so that, if one makes the replacement 
T a <S> T a — > 1(8)1, they become the corresponding color-singlet operators. The P-wave 
operators are normalized so that, if one makes the replacement (— |Pj)<8(— f Pj) — ► <%(1<8>1), 
they become (2j + l)/(2s + 1) times the corresponding S-wave operators. Given these 
normalization conventions, we find that the CEM relation (8) implies that the standard 
S'-wave and P-wave NRQCD matrix elements are related by 

(A*) = , 3(2j + 1) n n k 2 ^ x (o?Cs )}, (9) 

x n/ (2/ + l)(2/ + 3) maxX 1 v ;/ v ; 

where C n = 1 or (N 2 — 1)/(2N C ) = 4/3 if 0% is a color-singlet or color-octet operator, 
respectively. To be more specific, the CEM relation (8) implies that the S-wave matrix 
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elements are all related by simple group theory factors: 

(0f( s Si)> = 3(Of(^ )>, (10) 
{Of? S )) = ±(0*CSo)), (11) 
(Of ( 3 ^)> = 4<0f (%)>. (12) 

The P-wave matrix elements also include a factor /c^ ax from the integral over the relative 
momentum: 

(O^CPi)) = lk 2 max (0?eSo)), (13) 

(on 3 ^)) = ^*L*(o?Cs )), (i4) 
(ofCPi)) = l*L*(o?Cs )), (is) 

(Of? PA) = ^^ki^CSo)). (16) 

IV. ANALYSIS OF S- WAVE CHARMONIUM PRODUCTION 

In the production of 5-wave charmonium at the Tevatron with transverse momentum 
Pt > 5 GeV, it is known phenomenologically that the most important NRQCD matrix 
elements for H — J/ijj or ip(2S) are the color-octet matrix element (0%( 3 Si)) and a specific 
linear combination of color-octet matrix elements (Of ( 3 P )) and (O^^So)): 

Mf = (r/m 2 )((9f ( 3 P )> + (O?CS )), (17) 

where r m 3. Let us examine the ratio of these matrix elements 

off _ Ml (1X1 

where if stands for J/^ or ip(2S). The relation (9) yields the CEM ratio 

nH _ lvl r _ ' Sax _| L / 1Q \ 

^CEM" { QH {3Si)) ~ 15 m 2 + 3 - U»J 

The velocity-scaling rules of NRQCD predict that all of the matrix elements in both the 
numerator and denominator of Eq. (18) scale as v A relative to (0^( 3 Si)). Hence, the ratio 
in Eq. (18) is predicted to scale as v°. Since /c max scales as mv , the second term in the CEM 
ratio in Eq. (19) satisfies this scaling relation, but the first term does not. 



A. Analysis of Tevatron Data on J/ip Production 

We now use the relation (9) to test the validity of the CEM in comparisons with 
data. In Table II, we show values of R J ^ for several different sets of NRQCD ma- 
trix elements that have been extracted, under varying assumptions, from the transverse- 
momentum distribution of J/^'s produced at the Tevatron. The matrix elements were 
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TABLE II: Values of R J ^ , as denned in Eq. (18), in the NRQCD factorization approach and in the 
CEM. The column labeled u R J /1> n gives phenomenological values of R J ^ from various extractions 
of the NRQCD matrix elements from the CDF data [26]. The column labeled "Reference" gives the 
reference for each extraction, and the column labeled "PDF" gives the parton distribution that was 
used in the extraction. The headings "LO collinear factorization," "parton-shower radiation," and 
u kx smearing" refer to the method that was used to compute the NRQCD factorization prediction. 
The column labeled "-Rcem" gi yes the CEM ratios from Eq. (19) for the values of r and m c that 
were used in the NRQCD extractions of R J/lp . 



Reference 


PDF 


R J/* 


^CEM 


r 


m c (GeV) 


(k T ) (GeV) 


LO collinear factorization 


[27] 


MRS (DO) [28] 


10 ± 4 


0.44 


3 


1.48 






CTEQ4L [29] 


4.1 ± 1.2 tH 










[30] 


GRV-LO(94) [31] 
MRS(R2) [32] 


3.5 ± 1.1 ±1$ 

7.8 ± 1.9 ±fj;{j 


0.46 


3.5 


1.5 




[33] 


MRST-LO(98) [34] 
CTEQ5L [35] 


20 ± 4 
17 ± 4 


0.46 


3.4 


1.5 




parton-shower radiation 




CTEQ2L [36] 


1.4 ± 0.3 










[37] 


MRS (DO) [28] 
GRV-HO(94) [31] 


1.9 ± 0.6 
0.49 ± 0.11 


0.44 


3 


1.48 




[38] 


CTEQ4M [29] 


2.1 ± 0.8 


0.45 


3.5 


1.55 




kx smearing 


[39] 


CTEQ4M [29] 


5.7 ± 1.6 
2.6 ± 0.9 


0.46 


3.5 


1.5 


1.0 
1.5 


[40] 


MRS(D'_) [28] 


6.3 ± 1.7 
4.7 ± 1.2 


« 0.44 


3 


« 1.5 


0.7 
1.0 



taken from the compilation of Ref. [41]. We also show the CEM values -Rcem> taking 
fc^/m 2 = (m 2 D — m 2 c )/m 2 c ~ 0.54. Several sets of matrix elements were extracted by mak- 
ing use of NRQCD short-distance coefficients that were calculated at leading order in a s and 
under the assumption of standard collinear factorization. For these sets of matrix elements, 
R J I^ is much larger than Rq^ m - Multiple gluon radiation, as modeled by parton-shower 
Monte Carlos, tends to increase the partonic cross section more at smaller values of px than 
at larger values of px- Since the contribution of (0$^( 3 Si)) dominates that of M/^ at 
large px, while the contribution of ikf/^ is the more important one at small px, the effect 
of parton showering is to decrease the size of R J ^ . Hence, the addition of parton showering 
to the leading-order calculation of the NRQCD short- distance coefficients brings the ratio 
R J I^ for the extracted values of the matrix elements into better agreement with the -Rqem- 
However, the extraction that is based on the more recent CTEQ(4M) parton distributions is 
still in significant disagreement with -R^em- Surprisingly, kx smearing does not decrease the 
size of R J I^ as much as parton showering, and there is a substantial disagreement between 
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i?CEM an< ^ the v& l ues of R J ^ that are obtained by using kx smearing. 

Now let us compare the predictions of the CEM and NRQCD factorization directly with 
the CDF data [26]. 

The CEM predictions are from a calculation by Vogt [42] that makes use of the order-a;? 
cross section for production of a QQ pair [43]. Details of this calculation are given in Ref. [4]. 
The CEM factors Fr in Eq. (4) were fixed by comparison with fixed-target data. The charm- 
quark mass m c was also tuned to optimize the fits. The factorization and renormalization 
scales were chosen to be /i oc = \/ m c + p\, where pr is the sum of the transverse 
momenta of the Q and the Q. In our comparisons, we make use of the parameter sets labeled 
"■01" and "■04" in Ref. [4], which correspond, respectively, to fi — 2rriT, m c = 1.2 GeV, 
MRST98 HO [34] parton distributions and fi = m T , m c = 1.3 GeV, GRV98 HO [44] parton 
distributions. It is for these two parameter sets that CEM predictions of the pr distributions 
of charmonia produced at the Tevatron are available. 

The NRQCD predictions were generated from modified versions of computer codes cre- 
ated by Maltoni, Mangano, and Petrelli [45]. The codes compute the order-a 3 quarkonium 
production cross sections [46] and the standard Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution of the fragmentation contribution to the evolution of a QQ pair in a 
3 S\ color-octet state into a quarkonium. This fragmentation contribution is the dominant 
contribution at large pt- We chose m c = 1.5 GeV and took the factorization and renormal- 
ization scales to be /j, — rnr = \/ m c + Pt- O ur calculations made use of the same parton 
distributions as the CEM predictions, namely the MRST98 HO [34] parton distributions 
and the GRV98 HO [44] parton distributions. 

Plots of (data — theory) /theory are shown for the CEM and NRQCD factorization pre- 
dictions in Fig. 1. There is a substantial disagreement between the CEM predictions and 
the data. The normalizations of the CEM predictions are too small, and the slopes are 
relatively too positive. This discrepancy in the slopes is consistent with the fact that the 
CEM relation (19) over-estimates the size of {O^* ( 3 SJ) relative to . The NRQCD 
factorization predictions are in much better agreement with the data than the CEM predic- 
tions. Even if one were to adjust the normalizations of the CEM predictions to improve the 
fits, they would still be unsatisfactory, owing to the differences in slope between the CEM 
predictions and the data. 

A simple phenomenological model for the effects of multiple gluon emission on the the- 
oretical predictions is kx smearing. In kx smearing, the colliding partons are given Gaus- 
sian distributions in the intrinsic transverse momentum, with a width that is treated as 
a phenomenological parameter. A particular version of this model that has been used in 
comparing the CEM predictions with the experimental data [4, 42] attempts to account 
for multiple gluon emission from the two initial-state partons by adding two transverse- 
momentum "kicks" to the quarkonium momentum. The direction of each momentum kick is 
symmetrically distributed over the 4n solid angle, and the magnitude kx of each momentum 
kick is distributed as 

g(k T ) = ^e W (-k T /(k T )), (20) 

where (kj) is a phenomenological parameter. In the case of the CEM predictions with 
smearing, (/cf.) has been tuned to the value (fcf.) = 2.5 GeV 2 in order to obtain the best fit 
to the CDF J/ip data [4, 42]. This same value was also used in making CEM predictions for 
ip(2S) and Xc production. For purposes of comparison, we also take (kj) = 2.5 GeV 2 when 
we apply kr smearing to the NRQCD factorization predictions for charmonium production. 
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FIG. 1: J /ip production: (data — theory) /theory. The data are from the measurements of the 
CDF collaboration [26]. The upper figures are for the CEM predictions and the lower figures are 
for the NRQCD factorization predictions. The theoretical predictions for the left-hand and right- 
hand figures are based on the MRST98 HO [34] and the GRV98 HO [44] parton distributions, 
respectively. 

Plots of (data — theory) /theory for the CEM and NRQCD factorization predictions with 
kx smearing are shown in Fig. 2. The /^-smearing procedure substantially improves both 
the slope and normalization of the CEM fits to the data and slightly worsens the NRQCD 
factorization fits to the data. 3 The improvement of the slopes of the CEM predictions with 
kr smearing is to be expected since, generally, the effect of kr smearing is to increase the 



The effects of k T smearing on predictions for J/ip production cross sections at the Tevatron have also 
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FIG. 2: J/V' production: (data — theory) /theory, with fcr smearing of the theory, as described in 
the text. The data are from the measurements of the CDF collaboration [26]. The upper figures 
are for the CEM predictions and the lower figures are for the NRQCD factorization predictions. 
The theoretical predictions for the left-hand and right-hand figures are based on the MRST98 HO 
[34] and the GRV98 HO [44] parton distributions, respectively. 



been studied by Sridhar, Martin, and Stirling [40] and Pctrclli [39]. These studies made use of somewhat 
smaller values of (fcf.) than in the present work. They also concluded that the quality of the NRQCD 
factorization fits to the CDF data is little alfected by kr smearing. 
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TABLE III: Values of matrix elements, R J ^, and x 2 /d-o.f. from the NRQCD factorization and 
CEM fits to the J/ip data. In the NRQCD factorization fits, we set (Of /V, ( 3 5i))=1.16 GeV 3 and 
give the fitted values of (0^( 3 5i)) and . 



PDF 


<0 8 J/ *( s Si)> 
(GeV 3 x 10" 2 ) 


(GeV 3 x 10~ 2 ) 




X 2 /d.o.f. 


NRQCD Factorization 


MRST98 HO 


1.00 ± 0.22 


8.83 ± 1.24 


8.83 ± 2.27 


7.16/(ll-2)=0.80 


GRV98 HO 


1.02 ± 0.23 


10.6 ± 1.42 


10.4 ± 2.76 


7.98/(ll-2)=0.89 


MRST98 HO (smeared) 


1.41 ± 0.13 


0.41 ± 0.15 


0.29 ± 0.11 


10.28/(11-2)=1.14 


GRV98 HO (smeared) 


1.54 ± 0.14 


0.49 ± 0.16 


0.32 ± 0.11 


12.69/(11-2)=1.41 


Color-Evaporation Model 


MRST98 HO 








89.18/11=8.11 


GRV98 HO 








80.86/11=7.35 


MRST98 HO (smeared) 








20.78/(ll-l)=2.08 


GRV98 HO (smeared) 








45.70/(ll-l)=4.57 



cross section considerably at moderate values of pr and increase the cross section by a 
smaller amount at high p T . Nevertheless, the /cr-smeared CEM predictions yield poor fits 
to the data, as the slopes are still too positive relative to the data. This is consistent with 
the fact that the CEM relation (19) over-estimates the size of (0^( 3 Si)) relative to ikf/^, 
even for A/p-smeared extractions of the matrix elements. 

A compilation of values of matrix elements, values of R J ^, and chi-squared per degree 
of freedom (x 2 /d.o.f.) from the NRQCD factorization fits and CEM fits to the J/ip data is 
given in Table III. In the NRQCD factorization fits, the degrees of freedom are reduced by 
two, owing to the two NRQCD matrix elements that are varied in the fits. The unsmeared 
CEM fits have no free parameters, as the overall normalization is fixed by comparison with 
the fixed-target data. The smeared CEM fits to the J/ip data have one free parameter, 
which is then held constant in fits to the ip(2S) and Xc data. The values of R J ^ in 
Table III are much greater than Rq^ m ~ 0.46 in the fits without kx smearing and somewhat 

greater than -R^em m the n ^ s with kx smearing. The relative values of R J ^ and -Rqem are 
consistent with the large discrepancy between the slopes of the NRQCD factorization and 
CEM fits without k T smearing and the smaller discrepancy between the slopes of the NRQCD 
factorization and CEM fits with fc T smearing. 



B. Analysis of Tevatron Data on tp(2S) Production 

Next let us examine the case of ip(2S) production. In Table IV we show values of R^( 2S ) 
that were obtained from several different sets of NRQCD matrix elements that have been 
extracted from the transverse-momentum distribution of ip(2S) 's produced at the Tevatron. 
Again, the matrix elements were taken from the compilation of Ref. [41]. We also show the 
CEM values flggj in Table IV - As in the analysis of the J/ip data, the values of R^ 2S) lie 
substantially above Rq^m ■, except in the case of the matrix elements that were extracted 
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TABLE IV: Values of R^ 2S \ as defined in Eq. (18), in the NRQCD factorization approach and in 
the CEM. As in Table II, but for ip(2S). 



Reference 


PDF 


R f(2S) 


-"-CEM 


r 


m c (GeV) 


LO collinear factorization 


[27] 


MRS (DO) [28] 


3.8 ± 1.5 


0.44 


3 


1.48 


[30] 


CTEQ4L [29] 
GRV-LO(94) [31] 
MRS(R2) [32] 


4.1 ± 1.5 +H 
3.5 ± 1.3 tJJ 
7.8 ± 2.3 


0.46 


3.5 


1.5 


[33] 


MRST-LO(98) [34] 
CTEQ5L [35] 


3.1 ± 1.4 
2.1 ± 1.1 


0.46 


3.5 


1.5 


parton-shower radiation 


[37, 47] 


CTEQ2L [36] 
MRS (DO) [28] 
GRV-HO(94) [31] 


2.4 ± 0.8 

2.5 ± 0.9 
0.28 ± 0.35 


0.44 


3 


1.48 



by making use of the GRV-HO(94) parton distributions and parton-shower radiation in 
computing the NRQCD factorization predictions. 

Plots of (data — theory) /theory for the CEM and NRQCD factorization predictions with- 
out kx smearing for ip(2S) production are shown in Fig. 3. Owing to their larger error bars, 
the tp(2S) data have less discriminating power than the J/ip data. Nevertheless, it can be 
seen that the CEM predictions for ip(2S) production fit the data poorly, while the NRQCD 
factorization predictions fit the data almost too well. (It should be remembered that the 
error bars on the CDF data reflect both systematic and statistical uncertainties.) As in 
the case of J/ip production, the CEM predictions have overall normalizations that are too 
low and slopes that are too positive. Again, the discrepancies in the slopes are expected 
from the fact that the CEM relation (19) over-estimates the size of (Op 2S \ 3 Si)) relative to 
M$^ 2S \ Even if the normalizations of the CEM predictions were adjusted to fit the data, 
the fits would not be satisfactory. 

Plots of (data — theory) /theory for the /c^-smeared CEM and NRQCD factorization 
predictions for ip(2S) production are shown in Fig. 4. The effect of k T smearing is to 
improve both the CEM and NRQCD factorization fits. The CEM fits are improved in both 
normalization and slope and are now compatible with the data. The matrix elements and 
X 2 /d.o.f. from the NRQCD factorization fits and CEM fits to the ip(2S) data are given in 
Table V. The values of R^ 2S) in Table V are much greater than i?^ 2 ^ k, 0.46 in the 
fits without kr smearing and consistent with -Rcem* (and with zero) in the fits with kx 
smearing. The relative values of R^ 2S "> and -R^em are consistent with the large discrepancy 
between the slopes of the NRQCD factorization and CEM fits without kx smearing and the 
approximate agreement of the slopes of the NRQCD factorization and CEM fits with kx 
smearing. 
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FIG. 3: ip(2S) production: (data — theory) /theory. The plots are as in Fig. 1, but for ip(2S). 

V. ANALYSIS OF P-WAVE CHARMONIUM PRODUCTION 

Now let us turn to the case of production of the P-wave charmonium states Xcj {j = 0, 1, 2) 
at the Tevatron at transverse momenta pt > 5 GeV. It is known phenomenologically that the 
most important NRQCD matrix elements are the color-singlet matrix elements (Of C3 ( 3 Pj)) 
and the color-octet matrix elements (Og C3 ( 3 S'i)). The three color-singlet matrix elements 
can be expressed in terms of ((9^ c0 ( 3 Po)), and the three color-octet matrix elements can be 
expressed in terms of (Og c0 ( 3 S'i)) by making use of the heavy-quark spin-symmetry relations 

(^I( 3 ^)) = (2j + l)(^; 8 °( 3 Po)), (21) 
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FIG. 4: ip(2S) production: (data — theory) /theory, with kr smearing of the theory, as described 
in the text. The plots are as in Fig. 2, but for 4>{2S). 

which hold up to corrections of order v 2 . Therefore, we define a ratio 

((9f ( 3 Po))/m 2 ' 
The relation (9) yields the CEM prediction 



m 2 



R CEM - 1 5CfT5 • ( 23 ) 

The velocity-scaling rules of NRQCD predict that the ratio R Xc in Eq. (22) scales as v°. In 
contrast, we see that the CEM prediction in Eq. (23) scales as v ~ 2 . Furthermore, with the 
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TABLE V: Values of matrix elements, R^( 2S \ and x 2 /d-o.f. from the NRQCD factorization and 
CEM fits to the V(2S) data. In the NRQCD factorization fits, we set (Of (25) ( 3 5i))=0.76 GeV 3 
and give the fitted values of (of 2S) ( 3 Si)) and M^ S) . 



PDF 


{Ot {2S) ( 3 S 1 )) MtH S) R^ X 2 /d.o.f. 
(GeV 3 x 10~ 3 ) (GeV 3 x 10" 4 ) 


NRQCD Factorization 


MRST98 HO 
GRV98 HO 

MRST98 HO (smeared) 
GRV98 HO (smeared) 


2.34 ± 0.47 44.0 ± 19.2 18.83 ± 9.08 0.35/(ll-2)=0.04 
2.51 ± 0.52 55.4 ± 22.2 22.02 ± 9.93 0.55/(ll-2)=0.06 
2.12 ± 0.26 -6.77 ± 2.20 -3.19 ± 1.11 0.17/(ll-2)=0.02 
2.34 ± 0.29 -6.80 ± 2.39 -2.90 ± 1.08 0.22/(ll-2)=0.02 


Color-Evaporation Model 


MRST98 HO 
GRV98 HO 

MRST98 HO (smeared) 
GRV98 HO (smeared) 


47.72/11=4.34 
29.85/11=2.71 
10.43/11=0.95 
1.49/11=0.14 



TABLE VI: Values of R Xc , as defined in Eq. (22), in the NRQCD factorization approach and in 
the CEM. As in Table II, except for \c- The column labeled "-Rcem" § ives tne CEM ratios from 
Eq. (23) for the values of m c that were used in the NRQCD extractions of R Xc . 



Reference 


PDF 


R x - 


f>Xc 

-"-CEM 


m c (GeV) 


LO collinear factorization 


[27] 


MRS (DO) [28] 


(6.6 ± 0.8) xlO" 2 


36 


1.48 


[38] 


CTEQ4L [29] 


(0.71 ± 0.21) xl0~ 2 


40 


1.55 


[33] 


MRST-LO(98) [34] 
CTEQ5L [35] 


(5.8 ± l.l)xl0 -2 
(4.7 ± 0.8)xl0~ 2 


37 


1.5 



standard normalization of the NRQCD matrix elements in Ref. [1], the color factor in R Xc 
is estimated [46] to be 1/(2N C ) = 1/6, while the CEM prediction is that the color factor 
in R XQ is Cf = 4/3. Both the discrepancy in the velocity scaling and the discrepancy in 
the color factor have the effect of increasing the size of the CEM prediction relative to the 
expectation from NRQCD. 

In Table VI we show values of R Xc that were obtained from several different sets of 
NRQCD matrix elements that have been extracted from the transverse-momentum distri- 
bution of Xc's produced at the Tevatron. Again, the matrix elements were taken from the 
compilation of Ref. [41]. The values of R Xc in Table VI are reasonably close to the value 
R Xc ps v°/(2N c ) PS 0.17 that one would expect on the basis of the velocity-scaling rules and 
the estimate of the color factor. Table VI also contains the CEM values Rq c em - As expected, 
they are much larger than the values of R Xc that follow from the data. 4 Consequently, we 
expect the CEM to predict a cross section that is relatively too large at high p T , where the 



Because the CEM prediction scales incorrectly with v, we expect the disagreement between the CEM ratio 
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contribution from (Og c0 ( 3 Si)) dominates, in comparison with the cross section at low p?, 
where the contribution from (0 Xc0 ( 3 Po)) dominates. 

A compilation of values of matrix elements, values of R Xc , and x 2 /d.o.f. from the NRQCD 
factorization and CEM fits to the Xc data are given in Table VII. The values of R Xc in 
Table VII are much less than -Rqem ~ 37. 

Plots of (data — theory) /theory for the CEM and NRQCD factorization predictions with- 
out kx smearing for \c production are shown in Fig. 5. In the NRQCD factorization fits in the 
middle row, we have fixed the value of (O Xc0 ( 3 P )) to be 7.2 x 1(T 2 GeV 5 . This value is taken 
from a global fit to the existing data for inclusive decays of Xc states [48], which determines 
the corresponding decay matrix element to be (x c o|Ci( 3 -fo)|Xco) = (7.2 ± 0.9) x 1CT 2 GeV 5 . 
(Color-singlet decay and production matrix elements are simply related, up to corrections of 
order v 4 [1].) In the NRQCD factorization fits in the lower row, we have treated (O Xc0 ( 3 Pq)) 
as a free parameter. From the fits, it can be seen that both the CEM and 1-parameter 
NRQCD factorization predictions have normalizations that are too small and slopes that 
are too positive relative to the data. Both sets of fits show substantial disagreements with 
the slope of the data. However, the discrepancies are considerably greater in the CEM fits 
than in the NRQCD factorization fits. The large differences in slope between the NRQCD 
factorization and CEM predictions are consistent with the values of R Xc in Table VII rela- 
tive to -Rcem- The 2-parameter NRQCD factorization predictions fit the data much better 
than the 1-parameter predictions, but the values of (0 Xc0 ( 3 Po)) that are obtained are about 
a factor of six larger than the phenomenological value of the corresponding decay matrix 
element. Even given the large theoretical uncertainties in the determination of the decay 
matrix element, this discrepancy is unacceptably large. Thus, we conclude that the dis- 
agreement of the unsmeared NRQCD factorization prediction with the Xc data cannot be 
ameliorated by treating ((9f c0 ( 3 P )) as a free parameter. 

Plots of (data — theory) /theory for the /cy-smeared CEM and NRQCD factorization pre- 
dictions for Xc production are shown in Fig. 6. kr smearing improves the quality of the 
1-parameter NRQCD factorization fits and the quality of the CEM fit that is based on the 
MRST98 HO parton distributions, but does not substantially change the overall quality of 
the 2-parameter NRQCD factorization fits or the CEM fit that is based on the GRV98 HO 
parton distributions. The CEM predictions are still too large at high p T) and the 1-parameter 
NRQCD factorization predictions are now too large at low p^. The 2-parameter NRQCD 
factorization predictions with k? smearing give good fits to the data, but the fitted values 
of (Cj- c0 ( 3 Po)) now are considerably smaller than the central value from Xc decays. The 
relative values of R Xc and Rq% m are consistent with the relative slopes of the CEM and 
NRQCD predictions. The improvement in the 1-parameter NRQCD factorization fits with 
kr smearing suggests that multiple gluon emission may play an important role in quarko- 
nium production. However, the remaining discrepancy between the NRQCD factorization 
1-parameter fits and the data and the small value of (O Xc0 ( 3 P )) in the NRQCD factoriza- 
tion 2-parameter fits suggest that the value (kj.) = 2.5 GeV 2 is too large to be compatible 
with NRQCD factorization. In the fits to the J/ip and tp(2S) production data, one could 
compensate to some extent for changes in the value of (k^) by adjusting the value of M^ b . 
However, the corresponding adjustment of (0 Xc0 ( 3 Po)) in the fits to the Xc production data 



and the value extracted from the data to be even more dramatic for R Xb than for R Xc . For m b — 4.7 GeV 
and m B = 5.28 GeV, the CEM ratio predicted by Eq. (23) is Rq^ m ~ 76, while the expectation from 
NRQCD is that R Xb should be comparable to R Xc . 
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FIG. 5: Xc production: (data — theory) /theory. The plots are as in Fig. 1, except that they 
are for Xc the top row is for the CEM predictions, and the middle and bottom rows are for the 
1-parameter and 2-parameter NRQCD factorization predictions, respectively. 
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FIG. 6: Xc production: (data — theory) /theory, with kx smearing of the theory, as described in the 
text. The plots are as in Fig. 2, except that they are for Xc the top row is for the CEM predictions, 
and the middle and bottom rows are for the 1-parameter and 2-parameter NRQCD factorization 
predictions, respectively. 
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TABLE VII: Values of matrix elements, R Xc , and x 2 /d.o.f. from the NRQCD factorization and 



CEM fits to the Xc data. In the NRQCD factorization fits, the upper sets of parameters are for 
fits in which (0* cO ( 3 Po)) is fixed, as described in the text, while the lower sets of parameters are 
for fits in which (£>^ c0 ( 3 P )) is varied. 


PDF 


(0f°( 3 P o )) (Of°( 3 5i)) xVd-o.f. 
(GeV 5 x 1(T 2 ) (GeV 3 x 1(T 3 ) (1(T 2 ) 


NRQCD Factorization 


MRST98 HO 
GRV98 HO 

MRST98 HO (smeared) 
GRV98 HO (smeared) 


7.2 (input) 3.59 ± 0.39 11.23 ± 1.23 31.0/(11-1)=3.10 
7.2 (input) 3.94 ± 0.43 12.30 ± 1.35 35.5/(11- 1)=3.55 
7.2 (input) 1.71 ± 0.29 5.36 ± 0.89 17.4/(11- 1)=1.74 
7.2 (input) 2.08 ± 0.32 6.50 ± 0.99 14.5/(11- 1)=1.45 


MRST98 HO 
GRV98 HO 

MRST98 HO (smeared) 
GRV98 HO (smeared) 


40.8 ± 6.3 1.20 ± 0.60 0.66 ± 0.35 2.97/(ll-2)=0.33 
48.7 ± 7.3 1.17 ± 0.65 0.54 ± 0.31 3.19/(ll-2)=0.35 
3.88 ± 1.00 2.43 ± 0.36 14.12 ± 4.21 6.40/(ll-2)=0.71 
4.39 ± 1.09 2.67 ± 0.39 13.66 ± 3.93 7.88/(ll-2)=0.88 


Color-Evaporation Model 


MRST98 HO 
GRV98 HO 

MRST98 HO (smeared) 
GRV98 HO (smeared) 


50.20/11=4.56 
66.30/11=6.03 
16.15/11=1.47 
63.69/11=5.79 



is significantly constrained by the Xc decay data. Hence, comparisons of the predictions of 
NRQCD factorization with the Xc data may provide a more stringent test of NRQCD fac- 
torization than comparisons with the J/ip and ip{2S) data. A definitive test would require 
one to replace the /c T -smearing model for multiple gluon emission with a first-principles 
calculation. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have compared the CEM and NRQCD factorization approaches to in- 
clusive quarkonium production. As we have mentioned, the predictions of the CEM are at 
odds with a number of experimental observations. These include the different fractions of 
J/^'s from Xc decays that occur in B decays and in prompt production at the Tevatron, the 
nonzero polarization of J/^'s in e + e~ annihilation at the B factories, the nonzero polariza- 
tion of T(2S) and T(3S) in a fixed-target experiment, and the deviation from 3/5 of the 
ratio of the prompt-production cross sections for Xci and Xc2 at the Tevatron. Nevertheless, 
one might hope that the CEM would still be useful for predicting rates of inclusive quarko- 
nium production at large p T . While some of the predictions of the NRQCD factorization 
approach do not agree well with the data, for example, in the cases of the polarization of 
J/^'s produced at the Tevatron and J/ip photoproduction at the Hadron Electron Ring 
Accelerator (HERA) at Deutsches Elektronen- Synchrotron (DESY), the experimental and 
theoretical uncertainties are sufficiently large that one cannot yet make a definite statement 
about the validity of NRQCD factorization. (For a comprehensive review of these issues, 
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see Ref. [4].) 

By making use of the NRQCD factorization expressions for the quarkonium production 
cross section and for the perturbative QQ production cross section, we have translated the 
CEM assumptions into predictions for the ratios of NRQCD production matrix elements. 
In some instances, these predictions are at odds with the velocity-scaling rules of NRQCD. 
In other cases, they disagree with ratios of the nonperturbative NRQCD production matrix 
elements that have been extracted from phenomenology. Both of these facts indicate that 
the CEM picture for the evolution of a QQ pair into a quarkonium state is very different 
from that of NRQCD. 

A comparison of the CEM ratios with the phenomenological ratios that have been ex- 
tracted from the CDF data indicates that the CEM predicts a ratio M^ 1 / (0% ( 3 S\)) that 
is too small in J/ip and ip(2S) production and a ratio (Cg QO ( 3 1 Si))/(Ci go ( 3 P )) that is too 
large in Xc production. Both of these predictions of the CEM would be expected to lead to 
cross sections that have too positive a slope, as a function of pt, relative to the data. This 
expectation is borne out by comparisons of the CEM with the CDF data for J/ip, ip(2S), 
and Xc production. The CEM predictions do not yield satisfactory fits to the J/ip, ip{2S), 
or Xc data. The NRQCD factorization predictions yield satisfactory fits to the J/ip and 
ip(2S) data, but not to the Xc data, unless one relaxes the constraint on (O* ci) ( 3 P )} that 
follows from its relationship to the corresponding decay matrix element, which in turn is 
fixed through a global fit to the inclusive Xc decay data [48]. The normalizations of the 
CEM predictions are fixed through comparisons with the fixed-target data for charmonium 
production. 

kx smearing provides a phenomenological model for the effects of multiple gluon emission 
from the initial-state partons in a hard collision. Its effects are to smooth singularities at 
Pt = in fixed-order calculations, to increase the predicted cross section at moderately low 
Pt (away from the singular region), and to increase the predicted cross section by a smaller 
amount at high pt- Hence, the inclusion of A>r smearing would be expected to improve 
the fits of the CEM predictions to the charmonium data, which it does. Even with kT 
smearing, the CEM predictions show substantial disagreement with the data for J/ip and 
Xc production, but agree with the ip{2S) data, which have larger error bars. The smeared 
NRQCD factorization predictions are in good agreement with the data in the J/ip and ip(2S) 
cases and in reasonably good agreement in the Xc case. 

In general, we find that the nonperturbative NRQCD matrix elements can be adjusted so 
as to obtain good fits to the data, even when the slopes and normalizations of the partonic 
cross sections are modified quite strongly through the inclusion of kx smearing or the use 
of different parton distributions. Hence, the comparisons of the predictions of NRQCD 
factorization with the Tevatron data alone do not provide demanding tests of the theory. 
Comparisons with data from other processes are required. (See Ref. [4] for a review of such 
comparisons.) As we have mentioned, in the case of Xc production, the NRQCD factorization 
fits are constrained by the relationship of (0* cO ( 3 P o )) to the corresponding decay matrix 
element. Thus, there is less freedom in that case to tune the matrix elements to obtain 
a good fit to the data than in the cases of J/ip and ip(2S) production. Consequently, Xc 
production may provide a more stringent test of NRQCD factorization. The disagreement 
of the unsmeared NRQCD factorization prediction and the reasonable agreement of the 
smeared NRQCD factorization prediction with the shape of the Xc production data suggest 
that, if the NRQCD factorization picture is valid, then inclusion of the effects of multiple 
gluon emission is essential in obtaining the correct shape of the cross section. 



23 



Overall, the CEM predictions do not provide a satisfactory description of the data. Since 
the unsmeared CEM predictions for charmonium production at the Tevatron are absolute, 
one might reasonably argue that they could not be expected to fit the data as well as the 
NRQCD factorization predictions, which generally involve two free parameters. However, 
even if the normalizations of the CEM predictions are adjusted so as to optimize the fits, 
the predictions are still incompatible with the J/ip and Xc data, owing to incompatibilities 
between the predicted and observed slopes. In the case of the fc^-smeared predictions, the 
amount of kr smearing has been adjusted so as to optimize the fits to the J/ip data. Even 
so, the CEM predictions do not describe the data well. 

Finally, we should mention that there are several important sources of uncertainty in 
the theoretical predictions. These include uncertainties in the choices of factorization and 
renormalization scales, uncertainties in the parton distributions, uncertainties in the value of 
m c , and uncertainties from uncalculated corrections of higher order in a s . The uncertainty 
from uncalculated higher-order corrections is probably the largest. It is often estimated 
by comparing the differences between the predictions that are obtained by varying the 
factorization and renormalization scales by a factor of two. Such an exercise has been carried 
out in the case of the NRQCD factorization predictions for J/ip and ^>(2S) production [30] 
and suggests that the error from uncalculated higher-order corrections may be as large as 
100% at high pt and 40% at low pr- The CEM predictions presented in this paper involve 
a simultaneous change in scale and in parton distributions. If we assume that the change 
in scale is the dominant effect, then the uncertainty from it is about 75% at high p? and 
10% at low p T . It is likely that the large uncertainty in the NRQCD factorization prediction 
could be accommodated by a change in the values of the NRQCD matrix elements, except 
in the case of the Xc predictions. Even if we take the larger uncertainties in the NRQCD 
factorization prediction as being indicative of the uncertainties in the CEM predictions, 
that would not be enough to make the unsmeared CEM predictions compatible with the 
data. Therefore, it seems likely that, once the effects of multiple gluon emission are taken 
into account properly, through calculations of higher-order corrections rather than through 
the fcr-smearing model, the CEM predictions will still show a serious incompatibility with 
the data. As we have mentioned, calculations of the effects of multiple gluon emission are 
crucial to the NRQCD factorization predictions for Xc production, as well. Consequently, 
they could play an important role in sharpening the comparison between the CEM and 
NRQCD factorization. 
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